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We report on work studying the properties of the sphaleron in models of the elec- 
troweak interactions with two Higgs doublets in as model-independent a way as 
possible: by exploring the physical parameter space described by the masses and 
mixing angles of the Higgs particles. If one of the Higgs particles is heavy, there can 
be several sphaleron solutions, distinguished by their properties under parity and 
the behaviour of the Higgs field at the origin. In general, these solutions are not 
spherically symmetric, although the departure from spherical symmetry is small. 
SUSX-TH-01-001 

1 Introduction 

One of the major unsolved problems in particle cosmology is accounting for 
the baryon asymmetry of the Universe. This asymmetry is usually expressed 
in terms of the parameter 77, defined as the ratio between the baryon nunu 
ber density rig and the entropy density s: r\ — ub/s ~ 10 -10 . Sakharov u 
laid down the framework for any explanation: the theory of baryogenesis 
must contain B violation; C and CP violation; and a departure from_ther- 
mal equilibrium. All these conditions are met by the Standard Model □ and 
its extensions, and so there is considerable optimism that the origin of the 
baryon asymmetry can be found in physics accessible at current and planned 
accelerators (see a for reviews) . 

Current attention is focused on the Minimal Supersymmetric Standard 
Model, where_there are many sources of CP violation over and above the 
CKM matrix and the phase transition can be first order for Higgs masses 
up to 120 GeV, providing^the right-handed stop is very light and the left- 
handed stop very massive B. 

B violation is provided by sphalerons at a rate T s ~ exp(— E S (T)/T), 
where E S (T) is the energy of the sphaleron at temperature T. This rate must 
not be so large that the baryon asymmetry is removed behind the bubble wall, 
and this condition can be translated into a lower bound on the sphaleron mass 
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E S (T C )/T C <; 45. Thus it is clear that successful baryogenesis requires a careful 
calculation of the sphaleron properties. 

Here we report on work on sphalerons in the two-doublet Higgs model 
(2DHM) in which we study the properties of sphalerons in as general a set of 
realistic models as possible. In doing so we try to express parameter space in 
terms of physical quantities: Higgs masses and mixing angles, which helps us 
avoid regions of parameter space which have .akeady been ruled out by LEP. 

Previous work on sphalerons in 2DHMs llEralij has restricted either the 
Higgs potential or the ansatz in some way. Our potential is restricted only by a 
softly broken discrete symmetry imposed to minimize flavour-changing neutral 
currents (FCNCs). Our ansatz is the most general spherically symmetric one, 
including possible C and P violating field configurations EJ. 

We firstly check our results against the existing literature, principally 
BTT El, who found a new P-violating "relative winding" (RW) sphaleron, spe- 
cific to multi-doublet models, albeit at Ma =Mh± — 0. This is distinguished 
from Yaffe's P-violating deformed sphaleron or "bisphaleron" by a differ- 
ence in the behaviour of each of the two Higgs fields at the origin. We then 
reexamine the sphaleron in a more realistic part of parameter space, where 
Ma and Mh± are above their experimental bounds. We reiterate the point 
made in E3 that introducing Higgs sector CP violation makes a significant 
difference to the sphaleron mass (between ten and fifteen percent), and may 
significantly change bounds on the Higgs mass from electroweak baryogenesis. 



2 Two Higgs doublet electroweak theory 

The most general quartic potential for 2DHMs has 14 parameters, only one 
of which, the Higgs vacuuituexpectation value, v, is known. However, we are 
aided by the observation H that FCNCs can be suppressed by imposing a 
softly-broken discrete symmetry 0i(2) — > +(~)</ ) i(2)i an d results in a potential 
with 10 real parameters. One of these may be removed by a phase redefinition 
of the fields; the vacuum configuration is then entirely real, and CP violation 

is contained in one term 2^2 ^Re^l^) — ^j^J Im(</>j</) 2 ). Ignoring couplings 
to other fields, when \2 = there is a discrete symmetry (f> a — > — ier 2 0a, which 
sends Im^j^i-T - * — Im(</4</ , 2)- This can be identified as C invariance. 

Following E3 we determine as many as possible of the nine parameters 
in the potential from physical ones. The physical parameters at hand are 
the four masses of the Higgs particles, the three mixing angles of the neutral 
Higgses, one of which, 6c p, is the only CP violating physical parameter, (6c p 
mixes the CP even and CP odd neutral Higgs sector), and v. 
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We further check that the potential for these sets of physical parameters 
is always bounded from below. 



3 Sphaleron ansatz and numerical methods 

The most general static spherically symmetric ansatz is, in the radial gauge 



<t> a = (F a +iG a x a a a ) (^j 



r r 



•(i) 



Here, the subscript a = 1,2, and F a — a a + ib a and G a — c a + id a are 
complex functions. The boundary conditions can be most easily expressed in 
terms of the functions \i H a , h a , and Q a , defined by 

+ a a + ic a = H a e l0a , b a + id a = /i a e ie>Q , (2) 

and one can show that as r — > 0, either H% + l? a — > or Q\ — > ^>/2 + n\-K 
and 82 — > + n 2 7r, (ni,n 2 G Z). These boundary conditions distinguish 
between the various types of sphaleron solution: the ordinary sphaleron has 

+ Y? a — > as r — > 0, the bisphaleron has non- vanishing Higgs fields with 
ni = n-2, and the RW sphaleron non- vanishing Higgs fields with m 7^ n^- 
These integers represent the winding of the Higgs field around spheres of 
constant \4> a \i although only their difference has any gauge invariant meaning. 

Note that the ansatz is potentially inconsistent, as Im(<^>{<^2) oc £3, a point 
which has not been noted before. In practice, the non-spherically symmetric 
parts of the static energy functional, E[fA], contribute less than 1% of the 
total, and so we are justified in assuming the fields of the ansatz, /a, are a 
function only of the radial co-ordinate r and then integrating over £3 . 

We look for solutions to E[/a] using a Newton method which is an 
efficient way of finding extrema (and not just minima). The method can be 
briefly characterised as updating the fields Ja by 5 /a, given by the solution of 
E 5f = —E , where the primes denote functional diffcrentialion with respect 
to /a- A particular advantage to using this method is that because we calcu- 
late E , it is straightforward to get the curvature eigenvalues, and therefore 
to check that the solution really is the lowest energy unstable solution. 

4 Results 

We first checked our method and code against the results of BTT i, and 
Yaffetj finding agreement in the energy of better than 1 part in 10 3 for a wide 
range of parameters. We also measured the Chern-Simons numbers, ncs, of 
the solutions that they discovered and determined that they were near 1/2, 
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Energy of sphaleron and RWS 



Chern-Simons number of RWS 
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Figure 1. Contours of sphaleron energy (M\y/atw) an d Chern-Simons number. 



Most negative eigenvalue of sphaleron and RWS 



Second most negative eigenvalue of sphaleron 
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Figure 2. Eigenvalues (M m ) of the sphaleron solutions as a function of the CP even 
Higgs masses M h and Mjj. There was no mixing, and Ma = 241 Gev, M H ± = 161 Gev, 
tan/3 = 6, A3 = —0.05. The solid lines in the most negative eigenvalue plot represent the 
relative winding sphaleron, while the dashed lines are for the ordinary sphaleron, for the 
dotted region the potential was unbounded from below. 



but not exactly 1/2 as with the sphaleron solution. Further they appeared in 
P conjugate pairs with ncs of the pair adding to exactly one. We then looked 
at more realistic values of Ma and M H ± , with results that are displayed in 
Figs. |l] and [| Note first of all the well-known feature that the sphaleron mass 
depends mainly on Mh- Secondly, for increasing Mh, the curvature matrix of 
the sphaleron develops a second negative eigenvalue (Fig. ||, right), signalling 
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the appearance of a pair of RW sphalerons. The lower of the two ncs is 
plotted on the right Fig. [|. The departure from ncs = 1/2 is small, as is the 
difference in energy between the RW and ordinary sphalerons for the Higgs 
masses we examined, however the most negative curvature eigenvalue of the 
RW sphaleron can be double that of the sphaleron. More detailed results and 
discussion of their significance are reserved for a future publication LLa. 
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